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Radiation drag driven mass accretion in clumpy 
interstellar medium: implications for the supermassive 
black hole-to-bulge relation 
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ABSTRACT 

We quantitatively scrutinize the effects of the radiation drag arising from the radiation 
fields in a galactic bulge in order to examine the possibility that the radiation drag 
could be an effective mechanism to extract angular momentum in a spheroidal system 
like a bulge and allow plenty of gas to accrete onto the galactic center. For this purpose, 
we numerically solve the relativistic radiation hydrodynamical equation coupled with 
the accurate radiative transfer and quantitatively assess the radiation drag efficiency. 
As a result, we find that in an optically thick regime the radiation drag efficiency is 
sensitively dependent on the density distributions of interstellar medium (ISM). The 
efficiency drops according to in an optically thick uniform ISM, where tt is the 
total optical depth of the dusty ISM, whereas the efficiency remains almost constant 
at a high level if the ISM is clumpy. Hence, if the bulge formation begins with a star 
formation event in a clumpy ISM, the radiation drag will effectively work to remove 
the angular momentum and the accreted gas may form a supermassive black hole. As a 
natural consequence, this mechanism reproduces a putative linear relation between the 
mass of a supermassive black hole and the mass of a galactic bulge, although further 
detailed modeling for stellar evolution is required for the more precise prediction. 
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1 INTRODUCTION 

Recently, Kormendy & Richstone (1995) have pioneeringly 
suggested that the mass of a supermassive black hole (BH) 
does correlate linearly with the mass of the hosting bulge. 
(It is noted that the term of a bulge is used to mean a 
whole galaxy for an elliptical galaxy in this paper as is often 
so.) Further high-quality observations of the galactic center 
using stellar dynamics, gasdynamics, and maser dynamics 
(Miyoshi et al. 1995; Magorrian et al. 1998; Richstone et al. 
1998; Ho 1999; Wandel 1999; Kormendy & Ho 2000; Fer- 
rarese et al. 2000; Gebhardt et al. 2000a; Sarzi et al. 2001; 
Merritt & Ferrarese 2001a) allow us to make a detailed de- 
mography of supermassive BHs. The recent findings are the 
following: (1) The BH mass exhibits a linear relation to the 
bulge mass for a wide range of BH mass with a median BH 
mass fraction of /bh = A/ B H/M bu ig C = 0.001 - 0.006 (Kor- 
mendy & Richstone 1995; Richstone et al. 1998; Magorrian 
et al. 1998; Gebhardt et al. 2000a; Ferrarese & Merritt 2000; 
Merritt & Ferrarese 2001a). (2) The BH mass correlates with 
the velocity dispersion of bulge stars with a power-law re- 
lation as Mbh oc a n , n — 3.75 (Gebhardt et al. 2000a) or 
4.72 (Ferrarese & Merritt 2000; Merritt & Ferrarese 2001a, 



2001b). (3) /bh tends to grow with the age of youngest stars 
in a bulge until 10 9 yr (Merrifield, Forbes & Terlevich 2000). 
(4) In disc galaxies, the mass ratio is significantly smaller 
than 0.001 if the disc stars are included (Salucci et al. 2000; 
Sarzi et al. 2001). (5) For quasars, the /bh is on a simi- 
lar level to that for elliptical galaxies (Laor 1998; Mclure 
& Dunlop 2001a; Wandel 2001). (6) The /bh in Seyfert 1 
galaxies is under debate, which may be considerably smaller 
than 0.001 (Wandel 1999; Gebhardt et al. 2000a) or similar 
to that for ellipticals (McLure & Dunlop 2001a,2001b; Wan- 
del 2001), while the BH mass-to- velocity dispersion relation 
in Seyfert 1 galaxies seems to hold good in a similar way 
to elliptical galaxies (Gebhardt et al. 2000b; Nelson 2000; 
Ferrarese et al. 2001). These BH-to-bulge correlations sug- 
gest that the formation of a supermassive BH is physically 
connected with the formation of a galactic bulge. 

So far, very little is understood about the physical 
mechanism to produce such correlations, although some 
theoretical models has been proposed (Silk & Rees 1998; 
Ostriker 2000; Adams, Graff & Richstone 2001). Recently, 
as a possible mechanism to work in a spheroidal system, 
Umemura (2001) has considered the effects of radiation drag. 
The radiation drag is a relativistic effect, which may extract 
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the angular momentum effectively in a spheroidal system 
like a bulge, so that plenty of interstellar medium (ISM) 
could accrete on to the galactic center. Obviously, the radi- 
ation drag is inefficient in present-day elliptical galaxies or 
galactic bulges, since they possess little ISM. If the contents 
of a supermassive BH are initially in the form of ISM, how- 
ever, the bulge must have been optically thick in the early 
stage: 

where \ ls the mass extinction coefficient due to dust opac- 
ity, p is the density of the ISM, M gas is the mass of the ISM, 
and rb is the bulge radius. If a considerable amount of gas is 
expelled by a galactic wind at some stage, the optical depth 
should be still larger before the wind. If the radiation drag 
works efficiently in an optically thick medium, the rate of 
mass accretion induced by the radiation drag is maximally 
iboi/c 2 (Umemura, Fukue & Mineshige 1997, 1998; Fukue, 
Umemura & Mineshige 1997), where L b oi is the bolometric 
luminosity. Umemura (2001) has found that, if the maximal 
drag efficiency is achieved, the resultant BH-to-bulge mass 
ratio is basically determined by the energy conversion effi- 
ciency of the nuclear fusion from hydrogen to helium, i.e., 
0.007. However, it is not very clear whether this mechanism 
really works efficiently in realistic situations. 

In this paper, we investigate in detail the efficiency of 
the radiation drag in an optically thick ISM to test whether 
the radiation drag model is promising to account for the 
putative BH-to-bulge correlations. In particular, we concen- 
trate our attention on the effects of the inhomogeneity in 
the ISM. The model for the chemical evolution of elliptical 
galaxies suggests that an elliptical galaxy is initiated by a 
starburst in its early stage of (< 10 7 yr), and evolves pas- 
sively after a galactic wind event at a few x 10 s yr (Arimoto 
& Yoshii, 1986, 1987; Kodama & Arimoto, 1997; Mori et 
al. 1997). Also, in nearby starburst galaxies that have been 
studied, the ISM is highly clumpy (Sanders et al. 1988; Gor- 
don et al. 1997). Thus, if we consider the radiation drag in 
the early phase of bulge evolution, we should consider an 
inhomogeneous optically thick ISM. In this paper, to eluci- 
date the mutual effect between the dumpiness of ISM and 
the optical depth on the radiation drag efficiency, we build 
up a simple model of the bulge system and accurately solve 
the radiation transfer in a clumpy ISM. 

The paper is organized as follows. In Section 2, we con- 
struct the model of a galactic bulge. In Section 3, the basic 
equations for the ISM are provided. In Section 4, the angular 
momentum transfer efficiency is assessed in a uniform ISM. 
In Section 5, we investigate the angular momentum transfer 
efficiency by solving the radiation transfer in a clumpy ISM, 
and elucidate the relationship between the dumpiness of the 
ISM and the angular momentum transfer efficiency. In Sec- 
tion 6, we give implications for the correlation between the 
supermassive BH mass and the bulge mass. In addition, we 
discuss further effects that would give significant influence 
on the BH mass. Section 7 is devoted to our conclusions. 



2 MODEL 

We assume that a spherical galactic bulge consists of three 
components, that is, dark matter, stars, and a dusty ISM. 
The bulge radius r b is set to be l-10kpc by taking account 
of the observed sizes of elliptical galaxies or bulges in spi- 
ral galaxies. Dark matter is distributed uniformly inside the 
bulge. The mass of dark matter component, Mom, within 
the bulge size is equal to the stellar mass of the galactic 
bulge, Mbuige- As for the stellar component, we assume star 
clusters with a specific stellar initial mass function (see be- 
low). The star clusters are distributed uniformly inside the 
galactic bulge. Hereafter, 'a star' in this paper means 'a star 
cluster'. Here, AT*(= 100) stars are distributed randomly. 
For a dusty ISM, we consider the two cases: one is a uni- 
formly distributed ISM, and the other is a clumpy ISM. In 
the case of clumpy ISM, N c (— 10 4 ) identical clouds are dis- 
tributed randomly. The density p gas in a cloud is assumed 
to be uniform. The size of a gas cloud, r c , is a parame- 
ter. Then,the optical depth of a gas cloud is f = xPgas^c 
where \ is the mass extinction. We suppose that the stars 
and the ISM clouds corotate with the angular velocity cor- 
responding to the angular momentum obtained by the tidal 
torque at the linear stage of density fluctuations. Quantita- 
tively, the angular momentum is given by the spin parameter 
A = (Jt\E t \ 1/2 )/(GM t /2 ) = 0.05, where J T , E T , and M T 
are respectively the total angular momentum, energy, and 
mass (Barnes & Efstathiou 1987; Heavens & Peacock 1988). 
Here, the rigid rotation is assumed. 

The mass range of galactic bulges is postulated to be 
1O 6_13 M0. (However, in the present analysis, it is not very 
important to specify Mb u i ge , because the results are scaled 
with M bu i go as shown below.) The total mass of the ISM, 
M gas , and the mass of each gas cloud, m c , are parameters. If 
dust opacity as well as Thomson scattering is considered, the 
mass extinction is expressed by x = (icfT +n d a d ) / (p g + pa) , 
where o~t is the Thomson scattering cross section, n c is the 
electron number density, p g is the gas density, and rid, fa, 
and pd are respectively the number density, cross-section, 
and mass density of dust grains. If we take the dust-to-gas 
mass ratio /a g of ~ 10 -2 , then the opacity ratio is 
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where a<j is the grain radius and p s is the density of solid 
material within the grain. Thus, we find n^Cd n c UT in 
the situations of interest. Hence, in this paper, we evaluate 
the mass extinction by x — "d^d/ 'pg- 

Finally, as for the stellar evolution, we assume 
a Salpeter-type initial mass function (IMF) as cf> — 
A(m* I 'Mq) -1 for a mass range of [mi,m u ]. In the present 
analysis, we consider an initial starburst and the subsequent 
passive evolution of stars. The upper mass limit is inferred 
to be around 40M© in starburst regions (Doyon, Puxley & 
Joseph 1992). As for the lower mass limit, some authors 
suggest that the IMF in starburst regions is deficient in low- 
mass stars, with a cut-off of about 2 — 3M Q (Doane & Math- 
ews 1993; Chariot et al. 1993; Hill et al. 1994). In this paper, 
we assume m; = 2Mq and m u — 40Mq in the early stage 
of bulge formation. However, it should be also kept in mind 
that the lower mass limit is under debate; the claimed lower 
cut-off could be due to the magnitude limit effect (Selman et 
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al. 1999), and also recently sub-solar-mass stars are found 
in a starburst region in our galaxy, NGC 3603 (Brandl et 
al. 1999). In order to incorporate the stellar evolution, wc 
adopt the mass-luminosity relation (£,/Lq) — (m* / Mq) 3,7 , 
and the mass-age relation r* = 1.1 x lO lo (m„/M0)~ 2,7 yr 
(Lang 1974), where m„, £», and r* are respectively the stel- 
lar mass, luminosity, and age. 



calculate the radiation fields produced by spherically dis- 
tributed stars in an optically thin regime, and assess the 
angular momentum loss rate J and the mass accretion rate 
M. Next, we extend the analysis to an optically thick regime. 
Then, the relationship between the optical depth of a dusty 
ISM and the angular momentum transfer efficiency is de- 
rived. 



3 BASIC EQUATIONS 

As a relativistic result of radiative absorption and subse- 
quent re-emission, the radiation fields exert a drag force on 
moving material in resistance to its velocity. This radiation 
drag extracts angular momentum from the ISM, thereby al- 
lowing the gas to accrete on to the galactic center. The radi- 
ation drag is an effect of 0(v/c), but it could provide a key 
mechanism for angular momentum transfer in intense radia- 
tion fields. We put the origin at the center of the bulge, and 
adopt cylindrical coordinate r, <j), and z, where z-axis is the 
rotation axis of stars and gas. The components of the spe- 
cific radiation force that is exerted on moving fluid elements 
with velocity v are given by 



f T = ^(F r - v r E - Vl P" - v^P* 1 - « z P zr ), 



U = ^{F*- v^E - v+P++ - v r P r4, - v z P z4> ), 



and 

/.= 



-(F z - v z E - v z P zz - v T P IZ - Vrt> P 4 



(1) 

(2) 
(3) 



(Mihalas & Mihalas 1984) in the r-, 4>-, and z- directions 
respectively. Here, E is the radiation energy density, F a is 
the radiation flux, and P af3 is the radiation stress tensor 
where the non-diagonal components are null owing to the 
present symmetry. 

Using equations (1) — (3), we have the radiation hydro- 
dynamical equations to 0(v/c) as 



i*2S>_2=[F* -(E + P"*)v t ] 
r at c 



dt /g c [ 



(E + P zz )v z ] 



(4) 
(5) 
(6) 



where /g and /| are respectively r- and z- component of the 
gravitational force. 

The azimuthal equation of motion (5) is the equation of 
angular momentum transfer. This equation implies that the 
radiation flux force (the first term on the right-hand side) 
makes fluid elements tend to corotate with stars, whereas 
the radiation drag (the second term on the right-side) works 
to extract the angular momentum from gas. Therefore, the 
gain and loss of total angular momentum is determined by 
equation (5). 



4 UNIFORM ISM 

In this section, we consider the angular momentum transfer 
by radiation drag in a uniform ISM. First, we analytically 



4.1 Optically thin regime 

For uniform distributions of stars, the radiation fields inside 
the bulge are analytically integrated in an optically thin 
regime to be 
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where L hu \ gc and r b are the luminosity and radius of the 
bulge and lo is the angular velocity of stars. These quanti- 
ties are equivalent to those obtained by Fukue, Umemura 
& Mineshige (1997), except that the radiation flux F^ in 
the azimuthal directions is generated by the rotation of the 
bulge. As seen in equation (5), the flux F^ works to coro- 
tate the ISM with stars in contrast to the drag force. In Fig. 
1, we compare both forces exerted on gas per unit mass, 
jdrag = + pU^ v<j> / c and f = x F*/c, where = ruj. 
It is found that J drag overwhelms everywhere. Thus, the 
optically thin ISM can always lose the angular momentum 
as a result of radiation drag. The angular momentum loss 
rate per unit volume per unit time j is evaluated by 



3xw 



2tt 



■■f(x) 



(14) 



at the point of x = r/Vb, where f(x) — (^ + )x 2 + |a. 
Equation (14) is integrated over the volume of the bulge to 
give the total angular momentum loss rate J as 



J = 



•o-^bulge 3X^ 3 



(15) 



where po is the density of the ISM, and /3 = J xf(x)dV with 
dV being volume element in spherical coordinates. Noting 
that the initial angular momentum is 



2 2 
Jo = -M gas o)r b , 
o 

the mass accretion rate is expressed as 



sas Jo 
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(16) 



(17) 
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0.45 




Figure 1. The radiation drag force / dra s is compared with the 
flux force in the azimuthal directions. The i-axis is the radial 
distance normalized by the bulge radius r b , and the radiation 
force is normalized by (Sx^-t/bulgeV^TC 2 ?^). The figure shows 
that / dra s overwhelms everywhere. 



where r T = XPo r b = (3xA^gas)/(47rr b ). The coefficient r\ 
gives the radiation drag efficiency. By calculating the con- 
stant (3 numerically, rj is found to be 0.32. 



J _ n 2 (n + 4) 
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(23) 



where B = yfifT(n±l)/r(n±±), Ri is the mean distance from 
the rotational axis to each star, and T(n) is the Gamma 
function. It is noted that the difference between (22) and 
(23) is merely a small factor for ordinary rotation laws. As 
a result, it is found that the mass accretion rate decreases 
according to t^ 2 in an optically thick regime. The physical 
reason is that the larger is the optical depth of the bulge, the 
closer to zero is the difference between the velocity of stars 
and the gas velocity, so that the radiation drag efficiency 
falls. 

Combined with the optically thin case, we have the mass 
accretion rate in the uniform case as 



-^bulgc 



^bulgc 



(TT < 1), 

(tt > 1). 



(24) 



As readily understood by these results, the angular momen- 
tum transfer efficiency by the radiation drag is maximum 
when the optical depth of the ISM is around unity. 



5 CLUMPY ISM 



4.2 Optically thick regime 

In optically thick ISM, the regions where the optical depth r s 
from each star is less than unity are subject to the radiation 
drag (Tsuribe & Umemura 1997). In this situation, the flux 
for the region of r s < 1 is given by 
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(19) 



where R are R g are respectively the distance from the ro- 
tational axis to a star and to one volume element, dV = 
r 2 sin 9drd0d(f> in spherical coordinates for the position of a 
star. E = £ t /4ivcr 2 ,P 4 " t ' ~ with £„ being the luminosity of 
a star. Consequently, the local angular momentum transfer 
rate by equation (5) is 



£* 3 
cr 6 

This is summed up over all stars to give 

j _ 1 -^bulgc 



6 

where r s is the size of the region of r s = 1, and 
iV*£*. Then, the mass accretion rate is given as 
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(21) 
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So far, we have considered a rigidly rotating system. We can 
also analytically estimate the accretion rate in a system with 
a different rotation law as v$ ~ r n (e.g. n = —0.5 for the 
Keplarian rotation) . Then, the mass accretion rate turns out 
to be 



In the previous section, we considered the angular momen- 
tum transfer in a uniform ISM. In this section, we consider 
a clumpy ISM. First, we describe the treatment for the ex- 
tinction by clumpy gas clouds. Next, by solving the radiative 
transfer numerically, we assess the total angular momentum 
loss rate J and the mass accretion rate M gas . Then we elu- 
cidate the relation between the dumpiness of ISM and the 
angular momentum transfer efficiency quantitatively. 



5.1 Extinction by a clumpy ISM 

In the case of a clumpy ISM, N c (— 10 4 ) identical clouds 
with optical depth of f = xpgasfc are distributed randomly. 
In an optically thin regime, the radiation fields produced by 
a star are 



dE = 



c 47tr 2 



, dF„ 



47T7- 3 



r, dPo = dE , dP^ = dP§ z ~ 0,(25) 



where r = r^ — and r = |r| with i\ being the position of the 
i-th cloud and r, being the position of j-th star. The physical 
quantities with suffix are the radiation fields without ex- 
tinction. The corotating flux in the azimuthal direction and 
the flux contributing to the radiation drag are respectively 



= {dEo + dP$*)V„ 



dF£ ia * = -(dEo+dP$*)v* 



where V* = r,u) and v R 



(26) 
(27) 

s lj are the rotational velocities 



of a star and a gas cloud, with r* and r gas being respectively 
the distances from the rotational axis to a star and to a gas 
cloud. Next, we consider the extinction by dust in clumpy 
gas clouds. We calculate the radiation fields by the direct 
integration of the radiation transfer (see Fig. 2), where the 
extinction by all intervening clouds from a star to a cloud is 
summed up. Then, the radiation flux F, the radiation energy 
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density E, and the radiation stress tensor P are respectively 
given by 



F = ^dFojexp(-Tj), 
3=1 

iV* 

£ = ^dS , 3 exp(-Ti), 
3=1 

JV. 

P = y^dP ,jexp(-Tj), 



(28) 
(29) 
(30) 



3 = 1 



where tj = 2^)(1 - (6/r c ) 2 ) 1/2 r is the optical depth for all 
intervening clouds along the light ray. 



5.2 Angular momentum transfer in a clumpy ISM 

By substituting (26)- (30) for (5), we can evaluate the total 
angular momentum loss rate as 



J 



f; 



dragN. 



(31) 



where F" = ^ J x dF "j ex P(~ r j) with a being 'rot' or 
'drag', and n is the distance from the rotational axis to 
a gas cloud. If an optically thick cloud is irradiated by the 
radiation, only the optically thin surface layer is subject to 
radiation drag (Tsuribe & Umemura 1997). Then the sur- 
face layer is stripped by the drag and simultaneously loses 
angular momentum. We assume here that the stripped gas 
falls on to a central massive object. Then, with the above an- 
gular momentum loss rate, we estimate the total mass of the 
dusty ISM assembled on to a central massive object which 
may evolve into a supermassive BH. By the equation (31) 
and the relation M gas /M gas = — J/ Jo, the BH mass Mbh is 
assessed as 



rt(mi) _ rt(mi) 

= / M gas di = - / Mg, 
J o Jo 



'Jo 



dt, 



(32) 



where t(mi) is the age of a star with mass m;. 

Here, we introduce iVint as a measure of the dumpiness 
of ISM. The N int is defined by 



n c Tvr c r b = 



(33) 



where n c — N c /^irr b is the number density of gas clouds. 
iVint means the average number of gas clouds that are in- 
tersected by a light ray over a bulge radius. The degree of 
dumpiness is larger for smaller iVint, and 7Vint = oo corre- 
sponds to a uniform distribution. When 7V"i nt and the optical 
depth f of each cloud are specified, the total optical depth 
tt is given by tt = Ni nt f. With changing N- lnt , we inves- 
tigate the relationship between the dumpiness of ISM and 
the radiation drag efficiency. 

In Fig. 3, for different Mnt, we show the ratio of the 
resultant BH mass to the bulge mass (M B H/Mbui g e) against 
the total optical depth tt. In this figure, the dashed line 
is the analytic solution for a uniform ISM corresponding to 
JVint = oo, i e. equation (24). The solid lines show the nu- 
merical results for a clumpy ISM, where the arrows represent 
f = 1 for different Mnt • The thick lines show the results for 




Figure 2. A schematic illustration of the radiative transfer cal- 
culations in clumpy ISM. b and r c are respectively the impact 
parameter and the radius of the clumpy gas clouds, f is the opti- 
cal depth of a cloud. 



Miit > 1 and the thin lines are those for N- lnt < 1. In Fig. 3, 
several points are clear regarding the effects of the dumpi- 
ness on the radiation drag efficiency. First, if Mnt < 1, the 
drag efficiency is saturated when f > 1. Since the ISM is 
highly clumpy in this case, the radiation from distant sources 
which have different velocities from the absorbing clouds can 
contribute to the radiation drag even if tt 3> 1. However, the 
covering factor of clouds is smaller than unity, and therefore 
a large fraction of photons escape from the bulge without in- 
tersecting clouds. The saturation is attributed to these two 
effects. Secondly, if iVint ~ O(l), MsH/Afbuige grows with tt 
and maintains a high level in an optically thick ISM. Ob- 
viously, the behavior in the thick limit is different from the 
uniform case, although the situation is the same as the uni- 
form case in that almost all of photons are consumed inside 
the bulge. This is again because the dumpiness increases 
the mean free paths of photons, so that the radiation from 
distant sources enhances the drag efficiency. Thirdly, when 
TVint is much larger than unity, the distribution of ISM is 
closer to a uniform distributions and therefore the differ- 
ence between the velocity of a star and the velocity of an 
absorbing cloud is closer to zero, so that the drag efficiency 
falls as the uniform case. To conclude, the combined effects 
of the high dumpiness of the ISM and the large covering 
factor of clouds are necessary for high drag efficiency. 



6 BLACK HOLE-TO-BULGE MASS RATIO 

The linear relation between the BH mass and the bulge mass 
is a direct consequence of the present radiation hydrodynam- 
ical mechanism. The possible mass accreted by the radia- 
tion drag is given by M max = / ^ m '' Lbui S c(t)/c 2 dt. Incor- 
porating stellar evolution based on the present single star- 
burst model, we find M max = 1.2 x 10~ 3 Mb u i ge . As shown 
above, the drag efficiency is sensitively dependent on the 
density distributions of the ISM. In a highly clumpy ISM 
where JVint is smaller than a few, the BH-to-bulge mass ra- 
tio MBH/Afbui g o can be expressed by 



M m ax / 1 — tt \ 

Xf = V M 1 1 " 6 > 

Mbulge JWbulge 

= 1.2 x 10" 3 7?(1 - e" TT ) 



(34) 
(35) 



where rj gives the radiation drag efficiency. From the above 
analysis, rj is found to be maximally 0.34 in the optically 
thick limit, and then M B H/Mbui gc = 4.1 x 10~ 4 . 

The present model may be too simple in some respects 
to compare the results closely with the observations. For a 
more precise prediction of MBH/Mbui gc , it seems necessary 



6 N. Kawakatu and M. Umemura 




Figure 3. The BH-to-bulge mass ratio (MBH/M bulgo ) against 
the total optical depth (tt) of the bulge. The thick lines show 
the results for 7V; nt ^ 1 and the thin lines are those for JV; nt < 1. 
Filled circles denote JVi n t = 20, filled squares 7Vi nt = 5, filled 
triangles JV; nt = 1, open circles _/V; nt = 0.1, open squares iV; nt = 
0.01, and open triangles 7V; nt = 0.001. The dashed line is the 
analytic solution for a uniform ISM corresponding to JV; nt — > 
oo, where Mbh / M bu i ge oc tt in an optically thin regime and 
-^BH A^bulge oc TVjT in an optically thick regime. The arrows 
show the points where the optical depth of a cloud (f ) is unity. 
For JV; nt 1, the radiation drag efficiency is maximal if f 1. For 
Mnt < 1) the radiation drag efficiency is saturated when f ^ 1. 



to consider further effects which have not been incorporated 
in this simple model. In the present analysis, we assumed 
an initial coeval event of star formation, and the subsequent 
passive evolution without further star formation episodes. 
The radiation drag efficiency is basically determined by the 
total number of photons which are emitted from sources and 
absorbed by clouds during the whole history of the bulge. 
The recycling of the ISM for star formation generates more 
photons and therefore could enhance the mass ratio roughly 
by a factor of 2 (Umemura 2001). Also, in realistic situa- 
tions, the radiation drag is not likely to remove completely 
the angular momentum of stripped gas, and also stripped gas 
may be mixed with ISM having appreciable angular momen- 
tum. Although the detailed processes are not very clear at 
present, a little leftover angular momentum may lead to the 
formation of a viscous accretion disk around a BH, which 
could ignite QSO activity with nearly Eddington luminosity 
(Umemura 2001). If QSO activity is triggered, the radia- 
tion from the QSO can also exert the drag force. This effect 
can enhance the mass ratio by maximally a factor of 1.7 
(Umemura 2001). Furthermore, if a galactic wind occurs, it 
might cause a transition from an optically thick starburst 
phase to an optically thin QSO phase, and simultaneously 
reduce the drag efficiency. 

If the above effects are taken into account, the mass 
ratio could be maximally MBH/M bu i gc = 1.4 x 10~ 3 . This 
is compared with the recent observational data in Fig. 
4, combined with the result by a single-starburst model, 
MsH/Mbuigc = 4.1 x 10~ 4 . In Fig. 4, the hatched area is 
the prediction of the present analysis. As seen in this figure, 
the present radiation hydrodynamical model can roughly ac- 
count for the observational data, although the prediction 
falls slightly short of the observed median mass ratio. The 
prediction is smaller by a factor of six than (MeH/Mbuigc) = 
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Figure 4. The relation between the BH mass and the bulge 
mass. The vertical axis is the BH mass, and the horizontal axis is 
the bulge mass, in units of Mq. Recent observational results are 
plotted by symbols. The filled triangles denote elliptical galax- 
ies from Magorrian et al. (1998). The filled squares show ellip- 
tical galaxies from Ho (1999), Ferrarese & Merritt (2000), Ko- 
rmendy & Ho (2000), and Sarzi et al. (2001). The filled cir- 
cles denote elliptical galaxies from Merritt & Ferrarese (2001a). 
The small dots denote Seyfert galaxies from Ho (1999), Wandcl 
(1999), and Gebhardt et al. (2000a), and the open triangles show 
QSOs from Laor (1998). The relation of Magorrian et al. (1998) 
is Mbh = 0.006M bu[go , which is shown by a dot-dashed line; the 
relation of Merritt & Ferrarese (2001a) is M B h = 0.001M bulge , 
which is shown by a dashed line; and the relation for Seyfert 
galaxies is Mbh = 2.5 X 10~ 4 M bu i ge (Sarzi et al. 2001), which is 
shown by a thin dashed line. The hatched area shows the predic- 
tion of this papaer. The lower bound is a single-starburst model 
(Mbh = 4.1 X 10 — 4 M bu i gc ). The upper bound is the model in- 
coporating the effect of the recycling of star formation and AGN 
activity (M BH = 1.4 x 10- 3 M bulgo ). 



0.006 by Magorrian et al. (1998), while it is comparable to 
(M BH /M bu i gc ) = 0.001 by Merritt & Ferrarese (2001a). 

Another important effect is the geometrical dilution. 
In previous works, the radiation drag efficiency would be 
strongly subject to the effects of geometry (Umemura, 
Fukue, & Mineshige 1997, 1998; Ohsuga et al. 1999). A large 
fraction of emitted photons can escape from a disc-like sys- 
tem and thus the radiation drag efficiency is considerably re- 
duced. Recent observation have revealed that the BH mass 
fraction is significantly smaller than 0.001 in disc galaxies 
(Salucci et al. 2000; Sarzi et al. 2001). Geometrical dilution 
may be a reason for the observational fact in disc galax- 
ies, but the quantitative details are not clear before such an 
aspherical system is actually simulated. 



7 CONCLUSIONS 

By assuming a simple model of a bulge, we have investi- 
gated the mutual effect between the dumpiness of inter- 
stellar medium and the optical depth on the radiation drag 
efficiency for the angular momentum transfer. In a clumpy 
interstellar medium, we have accurately solved 3D radiation 
transfer to calculate the radiation drag force by the rotat- 
ing bulge stars. We find that the radiation drag efficiency 
is sensitively dependent on the density distribution of the 
ISM in an optically thick regime. The efficiency drops by a 
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factor of t,^ in a uniform ISM, while the efficiency turns 
out to be almost constant at a high level in a clumpy ISM. 
Also, the radiation drag efficiency falls as the covering factor 
of clouds becomes smaller than unity. Hence, for the radia- 
tion drag to work effectively, it is necessary that the cover- 
ing factor is larger than unity and that the distributions of 
ISM is highly clumpy. The present radiation hydrodynami- 
cal mechanism accounts for the linear relation between the 
mass of a supermassive black hole and the mass of a galac- 
tic bulge. The range of the predicted mass ratio is from 
M BH /Mbui g e = 4.1 x 10~ 4 to AfeH/A/buigc = 1.4 x 1(T 3 . 
Although we have adopted a highly simplified model, the 
result roughly accounts for the recent observational trends. 
Thus, it seems that more elaborate calculations deserve to 
be performed in the present context. In the future, we should 
consider further several effects which have not been incor- 
porated in the present analysis, and which can influence the 
mass ratio significantly: e.g. the realistic chemical evolution 
and geometrical dilution. 
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